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The density operator of the arbitrary physical system must be positive definite. Employing the 
general master equation technique which preserves this property we derive equations of motion for 
the density operator of an active atom which interacts collisionally with the reservoir of perturber 
atoms. The obtained general relations applied to the two-level atom yield Bloch-Boltzmann equa- 
tions (BBE). The form of the BBE obtained by us differs from that known from literature which, as 
we show, are not guaranteed to preserve the required positivity. We argue that our results are the 
correct ones and as such should be used in practical applications. Moreover, the structure and the 
terms which appear in our set of BBE seem to allow simpler and more straightforward physical 
interpretation. 
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I. INTRODUCTION 

A lot of experiments in atomic physics and spec- 
troscopy consists in investigating active atoms coupled 
to incident (laser) radiation and immersed in the ther- 
mal bath of perturbers which, typically, are the atoms of 
noble gas. Various phenomena occurring in such a sys- 
tem as well as its properties are then investigated. The 
amount of work devoted to such studies is enormous, it 
is therefore quite impossible even to list all relevant lit- 
erature, except a few essential monographs reviewing the 
subject |,§. 

The theoretical description of the discussed system 
must account for two major aspects. Firstly, the cou- 
pling of active atoms to the incoming radiation field 
and the radiative spontaneous phenomena must be prop- 
erly described. This is usually done by means of stan- 
dard methods of quantum optics [pi |j. As a result 
one obtains a set of equations of motion for the ma- 
trix elements of the atomic density operator. When the 
atom is described within a two-level model the obtained 
equations are known as optical Bloch equations. Sec- 
ondly, active atoms undergo collisions both with per- 
turbers and among themselves. The influence of colli- 
sions on physical properties of the system constituents is 
of paramount importance and is in itself a separate field 
of experimental and theoretical studies. Discussion of 
these problems in their full generality clearly goes beyond 
the scope of the present work. Let us however mention, 
that in spectroscopical applications the effect of collisions 
is usually accounted for by suitably derived (quantum- 
mechanically or classically) collision terms. Historically 
speaking, Boltzmann was the first to introduce the colli- 
sion terms into the equations of motion of the probability 
distributions. Therefore, for a two-level atom, the combi- 
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nation of optical Bloch equations together with collision 
terms might be called Bloch-Boltzmann equations (BBE) 
which account for both kinds of the discussed interactions 
influencing the behavior of active atoms. However, we 
feel it necessary to stress that the name Bloch-Boltzmann 
equations need not be restricted to two-level atoms. Gen- 
eralizations to more complex atomic models do not pose 
serious conceptual difficulties, though the form of corre- 
sponding equations of motion might be much more com- 
plicated. 

The main aim of this paper is to reexamine the origin 
and form of the collision terms in Bloch-Boltzmann equa- 
tions. We shall mostly study a simple two-level model, 
because it is formally the simplest, it allows the sim- 
plest interpretations, yet retaining the most important 
(at least in a qualitative manner) features of realistic 
physical situations and experiments. The motivation for 
our research is the following. First of all we note, that 
the density operator of an arbitrary physical system must 
always be positive definite. It is not clear whether the 
collision terms, used within the literature which is known 
to us, have this property. Moreover, they seem to exhibit 
other drawbacks or inconsistencies. These terms are de- 
rived (quantum-mcchanically or classically) by using ar- 
guments similar to those leading to collision integrals of 
standard Boltzmann equation pl| . We shall try to present 
a consistent theory which will, hopefully, allow us to clar- 
ify the question of positive definiteness as well as some 
more subtle points. 

The tools necessary to construct the proper form of 
BBE describing the system (active atoms) coupled to a 
reservoir (perturbers) are provided by the quantum the- 
ory of dynamical semigroups which entail the general 
master equation (ME) methods. It is worth stressing 
that we have in mind mathematically rigorous version of 
the ME theory based on completely positive quantum dy- 
namical semigroups. The extensive review of this subject 
is given in monographs [|5], where the authors derive 
and discuss the most general (sometimes called the Lind- 
blad form R pf) ME which preserves the positivity of the 
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considered density operator. Sec. II will, therefore, be 
devoted to a brief review of the essentials of this theory. 
In Sec. Ill we discuss how the general theory can be 
formally adapted to describe a mixture of two gaseous 
species when various models of internal structure of ac- 
tive atoms can be employed. We outline the procedures 
necessary to find explicit expressions for formal quan- 
tities introduced in the preceding section. We also try 
to identify physical assumptions needed to validate the 
presented theory. We restrict our attention to two-level 
atoms and, hence, in Sec. IV we give the derivation of 
the formal Bloch-Boltzmann equations in the above dis- 
cussed sense. Finally, Sec.V is devoted to the discussion 
of our results in the view of the facts known from the 
literature. Some features of the obtained BBE are dif- 
ferent from those known from other sources. We argue 
that our results are the proper ones. Since we employ 
much different theoretical formalism, the comparison of 
our results with other ones seems to be of an essential 
importance. 

As we have already stated, we focus our attention on 
the collisional phenomena occurring in the active-atom- 
perturber gaseous mixture. Therefore, we leave the ra- 
diative effects out of the picture. We shall, however, 
briefly indicate how, and under what conditions, such 
effects can be incorporated back into our formalism. 



II. THEORETICAL FRAMEWORK 

The physical system we will consider in this work, is 
a mixture of active atoms (A) with density Na and per- 
turbers of density N p . We assume that Na <C N p , which 
allows us to neglect the A — A collisions. Thus, only 
A — P collisions can affect the motion of the active atoms. 
Moreover, we assume that the density of perturbers is 
such, that only binary A — P collisions are of importance. 
The relatively dense perturber gas acts as a reservoir of 
energy and momentum and influences the velocity dis- 
tributions of the active atoms. On the other hand, it is 
usually safe to assume that P — P collisions are frequent 
enough to assure rapid thermalization of the perturbers. 
This allows us to think of a perturber bath as being in 
thermal equilibrium, so that the velocity distribution of 
P atoms is time-independent and given by a Maxwellian 

^W^^^exp^), (1) 

with = 2kBT/m p being the square of the most prob- 
able velocity of perturber atoms with mass m p , at tem- 
perature T. The stated physical conditions are not really 
very restrictive. They are fairly well satisfied by a great 
number of realistic experiments p| . 



A. Master equation for the quantum-classical 
system 

The master equation techniques stem from the theory 
of quantum dynamical semigroups applied to an open 
system, that is to a system which interacts with another 
one which serves as a reservoir. The active atoms may 
be considered as the quantum-mechanical open system 
which is coupled to the reservoir consisting of the per- 
turbers. The interaction between the two subsystems 
is manifested by the collisions occurring between A and 
P particles. The proper theoretical framework for the 
description of an open system, which ensures the preser- 
vation of positivity of the reduced density operator of the 
system of interest (active atoms), is supplied by the mas- 
ter equation (ME) approach j^]. It is not our aim to re- 
view the theory or derivation of ME. We shall rather ap- 
ply the general ME to a class of open systems which can 
be called quantum-classical ones. We study a quantum- 
mechanical open system, the states of which span the 
Hilbcrt space 

H = 0H a . (2) 

a 

The index a belongs to a set A which, for current pur- 
poses, is assumed to be discrete, but will subsequently 
be generalized to continuous one. A more physical in- 
terpretation of the spaces H a will be given later. We 
assume that within space Ti there exists a strong deco- 
herence mechanism [l^, |ll[ which practically excludes 
quantum superpositions of the form 

I ip ai ) ®\fp<x 2 ) ® (3) 

This is due to environmental decoherence mechanism, 
which is a generally accepted explanation of the ab- 
sence of superpositions of macroscopically distinguish- 
able states (Schrodinger cat problem) and emergence 
of classical properties, called dynamical superselection 
rules. The most effective physical mechanism leading to 
strong decoherence is provided by many subsequent colli- 
sions of the particles of the investigated system with the 
environment particles |^|. 

On each space H a we define a reduced density operator 
p a . The family of such operators 

P = {p a } a eA, (4) 

forms partially diagonal, quantum-classical density oper- 
ator which describes the properties of the relevant (open) 
system. The operators p a have the following properties 

(i) Pa ■ T~ta — > 'Ha] 

(ii) Pa > 0, positive — definiteness; 

(iii) V^Trp a = l, normalization. (5) 

a 

The most general form of the Markovian master 
equation for our quantum-classical density operator can 
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be obtained from the general form of the generator of 
the completely positive quantum dynamical semigroup 
(Lindblad-Gorini-Kossakowski-Sudarshan form). The 
relevance of complete positivity in the theory of quantum 
open systems is extensively discussed in Refs.[|[ ||. The 
corresponding master equation derived along these lines, 
which governs the evolution of the quantum-classical den- 
sity operator of the relevant system due to its interaction 
with the reservoir, reads (in the Schrodinger picture) 



dt 



£ 

B a Pa + PaB c 



(G) 



where additional index £ allows full flexibility to describe 
various dissipation phenomena. The operators intro- 
duced in this equation are defined as mappings: 

(i) H a = : H a — > H a , (hamiltonian) ; (7a) 

(ii) S^-.Hp^Ha; (7b) 

(iii) (»S'| a ) t : Ha -> Up, a map dual to S^; (7c) 

w ^ = EE(^) f ^- ( 7d ) 

The specific form of the Hamiltonian H a depends on the 
particular physical properties of the studied relevant sys- 
tem. Operators S^g depend on the interaction between 
the relevant subsystem and the reservoir. Recently the 
equation of this type has found its application in the 
quantum measurement theory (see the contribution of 
Blanchard and Jadczyk in where its properties are 
also widely discussed). 



B. Master equation for an atom immersed in 
perturber gas 

General and rather formal master equation (Q) has to 
be adapted to describe the presently discussed system, 
that is the moving (with velocity v) atom which collides 
with the perturbers. We shall proceed taking care of any 
additional or simplifying assumptions which go beyond 
the ones adopted in the derivation of ME (^). Moreover, 
we will proceed in a manner, which may be useful when 
considering active atoms with the internal structure more 
general than the simple two-level model. 

The collisions with perturber particles lead to strong 
decoherence which together with the uniform spatial dis- 
tribution of interacting particles justify the use of the 
density operators which are diagonal in momentum (or 
velocity) representation. Therefore, the discrete decom- 
position in Eq.(||) can be replaced by a continuous one 
which is taken to be with respect to the velocity v of an 
active atom. Thus, Eq.(B) is modified and becomes 



The ensemble of active particles is now described by a 
partially diagonal density operator p(v). We associate 
the space with the state space of an active atom which 
possesses velocity v. We introduce a set of operators 
{S^}, which constitutes a basis in the space of relevant 
operators acting on TL^. The specific form of operator ba- 
sis depends on the model chosen to describe the internal 
structure of an active atom. One may choose a multi- 
level model for which one has S a = Ski = k ) ( 1 1 , with 
k, I = 1,2,... , n, and with | k ) being the energy eigen- 
states. Alternatively, spherical tensor operators might be 
taken as a basis which is appropriate for atoms with spa- 
tially degenerate energy levels. Later on, we will consider 
a simple two-level model and we will explicitly define the 
necessary operator basis. 

First we analyze the hamiltonian term in ([|). To this 
end we expand it in the operator basis {S a }, and we write 
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H a — ► \h[w) =^> Q (v)Sa. 



(9) 



The particular form of the functions h a (y) need not be 
specified now. By the proper choice of these functions 
we can model various physical situations, some of which 
will be discussed later. For now, the first (hamiltonian) 
term of master equation becomes 



-i^2h a (v) 

a=l 



Sa, P(v) 



(10) 



It must be, however, noted that the Hamiltonian H(v) 
should be hermitian, so the functions h a (y) must satisfy 
some additional conditions, the particular form of which 
depend on the choice of the operator basis. We shall 
illustrate this point when applying the general formalism 
to derivation of the Bloch-Boltzmann equations for a two- 
level atom. 

Since the indices a and (3 are replaced by the "classi- 
cal" degrees of freedom, that is by velocities, when con- 
structing the second term of ME (||) we must replace the 
summation over the index (3 by integration over veloc- 
ities. Following the general rules given in 0] we may 
rewrite the second term in the ME as 

— E / Kab(?^v')S a p(v')Sl (11) 
a,b 

This term obviously has the sense of an operator which 
describes the transitions from a velocity group around 
v' to the velocity interval (v, v + dv). Hence it can be 
called a "gain" term. We shall postpone the discussion 
of the integral kernel to the further sections. At present, 
according to relation (|lTJ) , we shall only require that for 
any velocities v and v' 



U= dv H?, 



with 



v e 



(8) 



/C a h(v <— v ) — positively defined matrix, 



(12) 
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of the necessary dimensions. This matrix contains the 
details of the collisional interaction between the active 
and perturber atoms which will also be discussed later. 

Following further the principles of the construction of 
the ME m, we proceed to the third term in @. It is 
an anticommutator and it is built similarly to the former 
one. Namely, it can be rewritten as 



2 [j^aPa 



p a B a 



a , b 



v){StS b , p(v)}, (13) 



where the curly brackets denote an anticommutator. 
This term also describes the transition - escape from a 
velocity group (v, v + dv) to any other velocity so it is 
a "loss" term. 

Combining the discussed three terms, we now con- 
struct the master equation for a density operator of the 
moving active atom. We note, that no additional approx- 
imations (apart from those involved in the derivation of 
the general ME (ph) were made. Thus we have 



— p(v) = -i K{v) \s a , p(v) 



a, b 

■§W dv'/c: fc (v'^v){^5 b , P (m) 



a,b 



Let us note that the integration over velocity in the last 
term, affects only the integral kernel. Hence we can in- 
troduce the rate 



l* ab = 7 Q * 6 (v) = J dw 1 KTJ? «- v), 



(15) 



with the aid of which, our ME becomes 



^P(v)= - iJ2h a (v) S a , p(v) 



^Y,^){SlS b , p(v)} 



a . b 



J2 I 'dv' )C ab (v^v')S a p(v')Sl{W) 



The obtained ME is an operator equation. For practical 
purposes it is then convenient to expand the density op- 
erator p(v) in suitably selected operator basis. One can 
then compute all the necessary operator commutators 
and products, thus obtaining the equations of motion for 
the matrix elements of the density operator. We shall do 
so in the further section, by adopting a simple two-level 
model. Nevertheless, we stress that the presented ME 
( |l6| ) can be employed for atomic models more general 
than just the simple two-level one. 



III. MICROSCOPIC DERIVATION OF 

/C a6 (v^v') 

The formalism so far presented is fairly general. We 
proceed with its further discussion and clarification. The 
gain term ( [TT| ) and the loss one ( p"3| ) which describe the 
irreversible evolution of the relevant system stem from 
its interaction with environment. In our case, collisions 
are the manifestation of this interaction. The physical 
details concerning collisions are hidden in, so far rather 
formal, collision kernels K, ab {w <— v') which were left un- 
specified. Certainly their structure and mathematical 
properties follow from the procedures used when deriv- 
ing the necessary master equation. 

Derivation of the master equation for an open quantum 
system from the underlying fundamental Hamiltonian 
dynamics was the subject of very many investigations. 
Although the number of relevant literature sources is 
enormous, in only a few of them proper care is taken with 
respect to mathematical consistency of the presented re- 
sults. The density operator of an arbitrary system (in- 
teracting with the surroundings, or not) should be posi- 
tive definite. This can be ensured only by the carefully 
taken and properly conducted limiting procedures. It is 
not our aim to review these rigorous mathematical tech- 
niques such as weak coupling (or van Hove method) fll2[ , 
singular coupling |l3[ ] or low density limit fll4| , we refer 
the reader to [^) for a survey of the subject. 

We shall briefly discuss only the last of the mentioned 
limiting procedures - the low density limit which is de- 
signed specially for the description of a quantum system 
interacting collisionally with a low density perturber gas. 
The underlying physical assumptions are that only bi- 
nary collisions are of importance and that the duration 
of the collision is much shorter than the mean free-flight 
time (this latter condition is the essence of the so-called 
impact approximation). The reasoning leading to spe- 
cific form of the collision kernels is thus, as follows. Let 
us temporarily assume that the considered active atom 
is confined within a finite volume V and, therefore, is 
described by discrete states v, j ) with discrete veloci- 
ties and with quantum number j denoting its internal 
state (note that j may serve as a multiindex consisting 
of several quantum numbers). These states span certain 
Hilbert space, which by an introduction of an energy cut- 
off can be made finite dimensional. Let us say that N is 
the (finite) number of the basis vectors | v, j ). Then, we 
can use the completely positive generator obtained for 
a TV-state quantum system interacting with dilute gas 
of perturbers in the low density limit |l4| . Finally, one 
takes the density operator which is partially diagonal in 
velocities, removes the energy cut-off and takes the limit 
V — ► oo (a kind of a thermodynamic limit). Thus one 
arrives at the master equation which formally coincides 
with Eq. (|l6|) . The discussed procedure allows us to as- 
sign concrete meaning to all the terms which appear in 
the master equation. 



■5 



The Hamiltonian of the active atom is of the form 



the equation 



= Y d E j \j)(j\=Y,^j\j)u\ 



(17) 



where some eigenenergies Ej = Hujj may be degenerate 
and where n fixes the dimension of space of internal states 
of the active atom. This restricts the most general form 
of the Hamiltonian (J9J) to a diagonal one. Operators S a 
employed earlier, may be specified as eigensolutions to 



ff<°>, S a =hn a S a , a = 1,2, ,n 2 , (18) 



with £l a denoting the differences between atomic eigenfrc- 
quencies. Finally, the low density limit in the interaction 
terms leads to the following expression for the collision 
kernels 



^Cah(v 



(2tt) 4 h 2 



[LV r 



v v (v r - v r J 



dv' r I dv r S 3 



(19) 



where v (or v') are the velocities of an active atom, and 
v p (or Vp) of the perturber after (or before) collision. 
v r = v — v p (or v£ = v' — Vp) are the corresponding 
relative velocities, m, m p and /x denote the masses of 
the active atom, perturber and the reduced mass, re- 
spectively. Hp(v') describes the equilibrium velocity 
distribution of the perturber atoms, so in most cases 
it is just a Maxwellian as given in Eq.(Q). The func- 
tions T a (v, r <— v' r ) are related to the standard T-matrix 
(known from the quantum-mechanical scattering theory) 
in the center of mass variables 

n 

J2 T (v r ^v;)5 a = £<v r ,j|f \v' r ,f) \J)(J'\- 

(20) 



The resulting matrix /C a ;,(v <— v') is block-diagonal due 
to degeneracies f2 a = Qb which is ensured by the Kro- 
necker delta <5n a ,n 6 and obviously positively definite. 

For purposes of the further discussion it is convenient 
to express matrix /C a fc(v <— v') by standard scatter- 
ing amplitudes. From quantum scattering theory |L5| , 
we recall that the scattering amplitudes and functions 
T (v r <— vj.) are connected by the relation 

/»(v r - v' r ) = -(2tt) 2 A Tq(vV «_ v;). (21) 



Hence, rewriting Eq.([l9|) we alternatively have 



/Cab(v «- v') = ^ p <5o a ,0(, / dv^. I dv r 5 3 v - v' (v r - v' r 



x 5 



+ M7 a W p (v' - v' r ) f a (y r <- v' r ) / 6 *(v r <- v' r ). 



(22) 



Moreover, it might be also convenient to note that the authors. These forms are 
delta function reflecting energy conservation, can be writ- 
ten in several forms, each of them being used by different I ^v 2 




(23) 
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Equatio n (|l6|) , together with the Hamiltonian duced in we have 
specified in Eq.(|17]) and with collision kernels given by 
( |l9| ) or ( p2|) govern the evolution of the density operator 



of active atom interacting collisionally with perturbers. 
It preserves automatically positivity and normalization 
of the partially diagonal density operator p(v). 



IV. BLOCH-BOLTZMANN EQUATIONS 
A. Introductory remarks 

Our considerations presented in previous sections can 
be applied to multilevel atoms with complex internal 
structure. This is done by a suitable definition of op- 
erators S a introduced in Eq. (p8|) . Such a choice also 
specifies the form of the collision kernels. We shall, how- 
ever, illustrate the general approach by its application 
to a two-level atom which is a typical model for many 
quantum-optical phenomena. The resulting equations 
are, as we already mentioned, called Bloch-Boltzmann 
equations. In this case, the space TL^ of an atom having 
velocity v is isomorphic with C 2 . We denote the ground 
state by 1 1 ) and the excited state by | 2 ) . Taking them to 
form a basis in H^, we adopt the following identifications 



1) 



2) 



(24) 



The space of operators acting in is thus spanned by 
2x2 matrices. As a basis in the operator space we choose 
following four operators (pseudospin matrices) 



1) (1| = 

2) (2| = 

2><1 1 = 
1)(2| = 





1 

1 


1 




1 





-si 


(25a) 


= s 2 


- f?t 

— °2' 


(25b) 


= s + 


= 5*3 = $4) 


(25c) 


= S- 


= 5*4 = 5^ . 


(25d) 



The right-hand sides of these equations also fix the nota- 
tion we will use in this paper. For sake of easy reference 
we quote the nonvanishing commutators of operators S a : 



51, S3 

52, S3 



= -s 3 , 

= S3, 



51, S4 

52, S4 

53, S4 



= S4, 
= —S4, 
= 5*2 — Si. 



(26) 



Any operator on Ky can be expressed by the basis ones, 
while its matrix elements would be parameterized by ve- 
locity v. In particular, for a general Hamiltonian intro- 



~H(v) = hi(v)Si + h 2 ^)S 2 + h3(v)S3 + h i (v)S i 
n 



hi(v) fr 4 (v) 
h 3 (v) h 2 (y) 



(27) 



what is in accord with notation introduced in (Pq). We 
note, that the Hamiltonian must be Hermitian, which im- 
plies that functions hi(y), h 2 (v) are real, while /i3(v) = 
h%(v). Density operator p(v) can be expanded exactly in 
the same manner. To fix the notation and terminology 
we, therefore, write 

p(v) = Pll (v)|l)(l|+p 12 (v)|l)(2 

+p 21 (v)|2)(l|+p 22 (v)|2)(2| 
= pi(v)5i + p 2 (v)S 2 

+p 3 (v)^3 + i 04(v)^4. (28) 

Hence, we can identify the matrix elements 



Pi(v) = Pn(v), 

P3(v) =P2l(v), 



P2(v) = p22(v), (29a) 
p 4 (v) =p 12 (v), (29b) 



where the first row gives the populations, while the sec- 
ond one - coherences, which specify the physical meaning 
of the matrix elements of the considered density operator. 

Finally, we note that relation ([l7|) restricts the general 
form of the hamiltonian (|2^ ) to I13 = = 0, that is to 



: H (o) 



LOl 






(30) 



where Ej = hujj, (j = 1,2), are the eigenenergies of the 
corresponding levels of the active atom. Furthermore, 
it is straightforward to see that operators S a defined in 
( p5| ) satisfy relations JT§| ) with hamiltonian (pp|). The 
eigenvalue f^i = fl 2 = is doubly degenerate, while f^3 = 
lu 2 i, and f^4 = —lu 2 i (where Hu! 2 i = h(oj 2 — uji) > 0, is 
the energy difference between two atomic levels). 



B. General form of Bloch-Boltzmann equations 

Next step of our derivation consists in the expansion 
of the density operator of a two-level atom according to 
Eq. (p8|) . Doing so in the both sides of Eq.(|l6|) we then 
perform all the necessary operator computations. As a 
result we arrive at the following set of equations for each 
matrix element (|29| ) of the density operator 
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dt 



-ih 4 (v)p 3 (v) + ih 3 (v)p 4 (v) 

- (7l*l + 733) Pl(v) ~ \ [ill + 723) P3(v) - 5 (732 + 7*4) /°4(v) 



+ / rfv' 



P3(v) 



/Cn(v <- v')jOi(v') + /C 44 (v <- v')p 2 (v') 

+/C 4 i(v ^ v> 3 (v') + /C X4 (v <- v> 4 (v') 

i/i4(v)/9 3 (v) - ih 3 (v)pi(v) 

- (72*2 + 744) P2(v) - § (741 + 72*3) P3<v) - § (732 + 7*4) P4(v) 

+ y dw' [/C 33 (v v')Pi(v') + /C 22 (v v')P2(V) 

+/C 23 (v «- V> 3 (v') + ^32 (V <- v')P4(v') 

-ih 3 (y) (pi(v) - p 2 (v)) + i f/ii(v) - ft 2 (v)) p 3 (v) 
-3 (732 + 7*4) Pi W - 3 (^3*2 + 7i* 4 ) P2(V) 

-3 (7l*l + 722 + 733 + 74*4) P3(V) 



+ 



I" 



K 3 i{v <- v')pi(v') + /C 24 (v <- v> 2 (v') 

+/C 2 i(v <- v')/5 3 (v') + /C 34 (v <- v')p 4 (v') 



P4(v) 



i/i 4 (v) (pi(v) - /3a(v)j - i (/ii(v) - ?i2(v)J p 4 (v) 

-3 (72*3 + 7ll) Pl(v) - \ (72*3 + 74*l) P2(V) 

-| (7l*l + 72*2 + 733 + 7 44 ) P4(v) 



£i3(v <- v')jOi(v') + /C 42 (v <- v')p 2 (v') 

+/C 43 (v <- v')/03(v') + /Ci 2 (v <- v')p 4 (v') 

I 



(31a) 



(31b) 



(31c) 



(31d) 



This set of equations is the most general one, describing 
the evolution of active atoms (within a two-level model) 
due to the interaction with environment - collisions with 
the perturber atoms. Since these equations are derived 
directly from the general ME which preserves the pos- 
itivity of the reduced density operator of A atoms, we 
may be certain that this property is unchanged. 

The obtained BBE are somewhat simpler if we take 
into acco unt th e ha miltonian (|30|). The first lines of 
Eqs.(31aj) and (|31h| ) disappear, while the first lines of 
Eqs.( 31q ) and Q31d| ) become — iw 2 i/33(v) and iW2ip 4 (v), 
respectively. In the next subsections we will briefly 
present some frequently used approximations which al- 
low some simplifications of the general form (|3l]) with 
the first lines modified, as it was just discussed. We will 



first consider the so-called secular approximation and af- 
terwards we will assume that inelastic collisions can be 
neglected. However, it seems that the sequence of these 
approximations is irrelevant and that each of them can 
be used independently of the other one. 



C. Bloch-Boltzmann equations and secular 
approximation 

Secular approximation is recognized as a useful tool in 
the analysis of master equation technique applied to a 
manifold of physical systems. Its validity in quantum- 
optical problems is thoroughly discussed by Cohen- 
Tannoudji 0] and by Puri M. The essence of this ap- 
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proximation can, in the present context, be summarized 
as follows. When the off-diagonal terms in the Hamilto- 
nian (27) are much smaller (or simply negligible) as com- 
pared to the diagonal ones, the populations and the co- 
herences may be approximately decoupled. In the other 
words, under the secular approximation, the populations 
are coupled only to populations, while each of the coher- 



ences is coupled only to itself. So far, we have considered 
active atoms which are not subjected to any other in- 
teractions apart from the collisions. Since we have taken 
the hamiltonian (30), the secular approximation is clearly 
justified, which entails considerable simplification of the 
general Bloch-Boltzmann equations. Thus, we obtain 



dt 



Jt P ^ ] 



dt 



P4(v) 



- (7n(v)+73 3 (v))pi(v)+ / dv' [jC 11 (v^v')p 1 (v')+JC 4i (v^v')p 2 (v')\ , (32a) 

- (7 2 * 2 (v) + 74 * 4 (v)) p 2 (v) + J dv' [/C 33 (v <- v')pi(v') + /C 22 (v <- v> 2 (v')] , (32b) 
-icu 21 p 3 (v) - | ( 7 n(v) + 72 * 2 (v) + 73 * 3 (v) + 74*4 (v)) Pa(v) + J dv' £ 21 (v «- v> 3 (v"), (32c) 
iu} 2lPi {v) - \ (ju(v) + 7 22 (v) + 7 33 (v) + 7«(v)) pi(v) + / dv" /C 12 (v «- v')p 4 (v'). (32d) 



It should be noted, that if the active atoms are also irra- 
diated by incident radiation the Hamiltonian ( |30|) will be 
modified by the suitable coupling terms. Moreover, one 
has also to account for spontaneous emission, i.e., for the 
coupling to the vacuum field. This would lead to the ap- 
pearance of additional terms describing radiative effects. 
In such a case the validity of the secular approximation 
must be separately investigated. 



D. No inelastic collisions 

The energy transfer during the typical collision be- 
tween the atoms in the gaseous mixture is of the order of 
ksT. Since the typical temperatures of spectroscopical 

I 



r 



experiments are of the order of several hundred kelvins, 
the energy available during the collision is by two orders 
of magnitude smaller than the typical separation of the 
levels of the active atom. Hence, the probability that 
the collision would excite or deexcite the active atom is 
negligible. This means that the inelastic collisions can 
be left out of the picture. Inspection of Eqs.(|3l|) or ( |3^ ) 
tells us that relations 



for a = 3, 4 : K aa {v <- v') = 0, 



T* = 



(33) 



are equivalent to neglecting the inelastic collisions. Em- 
ploying such an assumption in Eqs.(J33) we get 



dt 



P2(v) 



dt 



P4(v) 



7ii(v>i(v)+ / dv' /C n (v<-v')pi(v'), 



-7 22 (v> 2 (v) + J dv' /C 22 (v<- v')p 2 (v'), 

-iw 2 ip 3 (v) - | (7n(v) +7 22 (v)) p 3 (v) + J dv' JC 2 i(v <- v')p 3 (v'), 
-iw 2 ip 4 (v) - \ (711 (v) + 722 (v)) p 4 (v) + / dv' /C i2 (v <- v')p4(v')- 



(34a) 
(34b) 
(34c) 
(34d) 



Bloch-Boltzmann equations in this form seem to be used Our derivation clarifies the procedures leading to these 
in quantum-optical applications which are known to us. 
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equations. It also specifies the methods allowing explicit 
computation both of collision kernels and rates. More- 
over, general form (|3l]) of the BBE ensures preservation 
of the positivity of the density operator of active atoms 
due to the conditions ( |l2| ) and ( |Tq ) which seem not to be 
used in the literature. 



V. DISCUSSION 
A. Review of the results known from literature 



cerned with a simple two-level atom suffering elastic colli- 
sions, we will restrict our attention only to such a model. 
The derivation of the collisional terms, as discussed in 
|l| is performed by means of a truncation of the BBGKY 
hierarchy of equations of motion to one- and two-particle 
density operators. Assuming that the active atoms and 
perturbers are uncorrelated for time t — > — oo, and treat- 
ing their translational motion quasi classically, one then 
arrives at the kinetic equation which governs the evolu- 
tion of the matrix elements of the density operator due 
to collisions. This equation is of the general form 



The theoretical models employed to describe atomic 
A — P collisions which tend to thermalize the velocity 
states of active atoms are usually based on the suit- 
able adaptation of the quantum-mechanical, or classical 
Boltzmann equation. The most extensive and thorough 
review of such an approach is given in a monograph by 
Rautian and Shalagin (!]]. These authors present quite 
general formalism which can be applied to describe var- 
ious physical situations concerning active atoms with 
complex, multilevel structure of degenerate energy eigen- 
states. Moreover, general expressions can also treat both 
elastic and inelastic collisions. Since we are mainly con- 



dt 



~ 7ii(v)Py(v) 

+ / dv' W y (v^v')^(v'),(35) 



for i,j = 1, 2. It is valid for binary elastic collisions in 
the impact approximation and for perturbers with ho- 
mogeneous spatial distribution. The collisional kernel 
Wjj(v <— v') following from quantum-mechanical Boltz- 
mann equation is expressed via the corresponding scat- 
tering amplitudes (see [Q, Eq.(2.147)) 



Wij(v <- v') = 2N p dv' r / dv r S 



— ->l r* 1 7* —I • 

v- v (v r -v r 

m 



6 [yl - tv 2 ) W p (-v' - w' r ) fi(v r <- w' r ) f*(v r <- v^.). 



(36) 



The notation employed here is exactly the same as that 
in Eq.([l9|), while /,(v r <— v'.) is the elastic scattering 
amplitude of an atom in the internal state | j_V The col- 
lision rate 7ij(v), known from literature fj], W i has the 
form 

7y (v) = N p ^— J dv r Wp(v - vy) 

x [/i(v r ^v r )-/;(v r ^v r )], (37) 

so it is given by the difference between the scattering 
amplitudes corresponding to forward scattering. 

The diagonal kernels Wu(v <— v') and rates ju are 
seen to be real and can be thought of to possess classical 
limits. Eq.(|35|) multiplied by dv (for i — j) represents 
the collisional diffusion out of and the flow into veloc- 
ity interval (v, v + dv) due to scattering of active atoms 
in \i) state. Therefore, integral term in Eq.(35) is the 
gain one and it gives the number of active atoms in state 
| i ) which change velocity from v' before, to v after the 
collision. Hence, the collision kernels Wa{v <— v') are 
measures of transition probabilities from v' to v velocity 
groups. The rates 7ii(v) enter the loss terms, and give 
the number of atoms in the state | i ) escaping from ve- 



locity interval (v, v + dv) to any other one. Let us note 
that 7ii(v) can be also viewed as the collision frequency 
and its inverse t^(v) = l/ju (v) can be interpreted as 
the average time between collisions. Finally, it is worth 
noting that the given probabilistic interpretation of the 
kernel and frequency is fully consistent with the relation 



7«(v) = / dvi W u (v<- v'), 



(38) 



which reflects the requirement of the particle number 
conservation and follows from Eqs.(|36|) and ( |37j ) by the 
application of the optical theorem. 

The quantum-mechanical (linear) Boltzmann equation 
( |35| ) specifies the time evolution not only of the diagonal 
elements of the density operator, but also of the coher- 
ences - the off-diagonal ones, that is for i ^ j. In such 
a case, neither the kernels Wy(v <— v') nor the rates jij 
have probabilistic interpretation. The kinetic equation 
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( p5| ) for coherences can be recast into the following form 



dt 



+ [<& Wijtv*- v')Pij(v')- (39) 

The collisional rate is split into two terms: 7y (v) 

being called the "phase-changing" rate, while 7]^ (v) 
is the "velocity changing" one Jl6|. The former rate, 

lij h \v), is the on c usually associated with the homo- 
geneous linewidth which appears in the typical pressure 
broadening theories of the atomic spectral lineshape. It 
can be argued |l6| ] that, in many practical applications, 
7iJ may be treated as velocity independent (although 
this dependence entails interesting correlations between 
Doppler and power broadening). Two last terms in (|39| ) 
describe velocity changing collisions and are frequently 
omitted from the equations of motion of density oper- 
ator. There are some arguments |l(| justifying such an 
omission. However, these arguments are somewhat vague 
and heuristic, and according to our knowledge, there is 
no clear-cut and generally accepted criterion which speci- 
fies when one is allowed to drop out the velocity changing 
terms in 

In many theoretical approaches used to describe the 
spectroscopical experiments the semiclassical model is 
thought to be sufficient. In such a case, the motion of 
atoms is taken to be purely classical, while the internal 
degrees of freedom are treated quantum-mechanically. 
Then, the classical kinetic theory is quite appropriate for 
determination of velocity distributions of not too dense 
gases. The discussed gas is then considered as a mixture 
of the following species: perturbers, active atoms in the 
ground state and atoms in the excited state. Since the 
perturbers are assumed to be in thermal equilibrium, one 
can invoke linear Boltzmann equation jl7j and derive ki- 
netic equation for classical velocity distributions Fj(v) 
which are (up to a normalization constant) proportional 
to the probability of an active atom being in state | j ) 
and possessing velocity v, that is, to populations pu(v). 
The resulting kinetic equation has the form identical to 
diagonal equations ((3^) . The classical kernel for an active 
atom in state \ j ) is then 



S vc h 



col I 



Wf'^v^) = N p I dv. 



5{v r - v' r ) 



v-v'--^(v r -v;)l, (40) 



and the corresponding classical collisional rate is 



7 ^(v) = N p I dv p / dfi(0) |v r | Wp(vp) 



doj(v r , d) 



dfl 



(41) 



These quantities also satisfy requirement ( p8| ) and have 
the discussed probabilistic interpretation. The collisional 



(classical, center-of-mass) cross sections daj/dft may ob- 
viously be viewed as classical counterparts of scattering 
amplitudes appearing in "diagonal" definitions ( ^6|) and 
(p7|). There are no classical analogs of coherences, so 
there is no classical argument leading to coherence ker- 
nels and rates. Atomic coherences describe the correla- 
tions between the internal quantum states of active atoms 
and usually may be left out of the classical picture. 

Kinetic equations ( |35| ) or their semiclassical counter- 
parts, is then used as an augmentations of the equations 
of motion describing the radiative effects occurring due 
to the coupling of the active atoms with the incident 
(laser) radiation. Such an approach is frequently used in 
the studies light-induced kinetic effects in gases [|l8[ |l^| 
and of other phenomena of the spectroscopic nature in 
which atom-atom collisions are of importance (e.g. p(i|). 
However, computation of the collision kernels and rates 
either in quantum-mechanical or in classical case is usu- 
ally extremely difficult, if at all possible. Calculations 
of the scattering amplitudes (or equivalently of the cross 
sections) require the knowledge of interatomic potentials. 
Such potentials are either not known, or given only in the 
numerical form. Therefore, many authors employ model 
kernels and rates which allow analytical computations. 
The most popularly used models are strong collision one 
and the Keilson-Storer model. There are also other pos- 
sibilities connected with the notion of collision operators. 
These models and concepts are reviewed elsewhere, see 
for example ]l|, |2l], |2^] and references given therein. 



B. Comparison to our results 

We are now in position to compare our results with 
those above presented which are known from literature. 
Although we focus our attention on the Bloch-Boltzmann 
equations, i.e., on those for a two-level atom, it seems 
that the generalizations to more complex atomic struc- 
tures are rather straightforward. 

Adopting the operator basis (^||) with the hamiltonian 
( p7| ) and leaving the inelastic collisions out of the picture, 
which corresponds to putting T3(vy <— v£.) = T±(v r «— 
v',) = 0, we arrive at the BBE exactly as in Eqs.(|34"|). 
Since we have an identification ( |29| ) we can conclude that 
our collision kernels /C a h(v <— v') coincide with Wij(v <— 
v'), for a — i and b — 7, they are given by the same 
formulas, Eqs.(|l9|) or ( f22|) and (|36|) , respectively. Hence, 
the comments given after Eq.(p7|) apply to population 
kernels Wu(v <— v') as well as to JC aa {^ <— v') for i — 
o = l,2, 

The situation is, however, much different with respect 
to our collisional rates 7* b as compared to literature ones 
jij. First of all, we note that in our case the rates 7* h 
are defined as integrals over kernels /C a (,(v <— v') and fol- 
low quite naturally from the general formalism. Hence, 
"fab are gi ven as integrals over the products of scattering 
amplitudes. The rates 7^ from literature are given by 
differences of forward scattering amplitudes, see Eq.([37|). 
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The derivation and the limiting procedures, as presented 
by Rautian and Shalagin, which lead to the expression 
( |37| ) seem to lack mathematical rigor and their physical 
basis seems to us to be not fully clear. Luckily, the appli- 
cation of the optical theorem entails correct relation (pq) 
which, on the other hand, is an automatic consequence 
of our formalism. 

Secondly, as it follows from the very structure of our 
theory, Bloch-Boltzmann equations for coherences, that 
is Eqs. (|34c[) and (34d) contain the coherence collisional 



rate, which is in our case, of the form 

ri2(v) = r 21 (V) = i[ 7ll (v) + 72*2(v 



(42) 



Such a quantity does not appear in the kinetic equation 
approach. The presence of Ti2 is connected with the em- 
ployed method which ensures the positivity of the atomic 
density operator. The kinetic equation method which 
provides coherence rate 7i 2 (v) cannot, therefore, guaran- 
tee that the computed density operator will be positively 
defined, as it always should be. 

Moreover, the rate r 12 (v) is specified only by col- 
lision kernels /Cn(v «— v') and /C 22 (v <— v') which 
may be called "population" ones. Fi 2 (v) is, thus, in- 
dependent of the coherence kernel /Ci 2 (v <— v') which 
enters Eqs.([}4|) separately. Only if inelastic collisions 
are taken into account, that is if the scattering ampli- 
tudes /s(v r <— v' r ) and fi{y r <— v' r ) are nonzero, the 
off-diagonal rates 7* h may appear in Bloch-Boltzmann 
equations. This is clearly seen by inspection of Eqs.(pl|) 
which are the most general ones. We are of the opin- 
ion that our Ti 2 (v) should be associated with homo- 
geneous linewidth, usually connected with the so-called 
"phase-changing" collisions. Having in view the relation 
between Ti 2 and population collision kernels, we sup- 
port the opinion of Rautian and Shalagin (jjj, p. 66), that 
any attempts to distinguish "phase-changing collisions", 
"velocity-changing collisions" or any other special types 
of collisions, are futile and physically misleading. Thus, 



it also seems that attempts to split the collisional rate as 
in Eq.(|39|) are arbitrary and not really justified. 

Furthermore, if we assume that the coherence kernel 
K%2 (v <— v') can be neglected (putting aside the ques- 
tion of justification of such an approximation) the rate 
r 12 (v) appears naturally in our formalism. There is no 
need to invoke any additional arguments that homoge- 
neous linewidth should be present in Bloch-Boltzmann 
equations, even when the kernel /Ci 2 (v <— v') is absent. 
This conclusion remains valid, also in the case we adopt 
semiclassical picture taking classical expressions ([ll]) and 
©for 7 * Q and/C oa (v^v'). 

Finally, we note that relation (|3^) which is the require- 
ment of particle number conservation follows automati- 
cally from our derivation, as it is reflected by expression 
( |l5| ) . This is especially important in various approximate 
approaches to estimate the collision kernels. In such ap- 
proximations the optical theorem may not hold, and thus 
( |38| ) may also not hold. The kinetic equation method 
may not ensure the particle number conservation while 
our approach clearly does, even if some approximation 
are made. 

We conclude this work, by saying that we believe that 
the proper form of the Bloch-Boltzmann equations is ex- 
pressed by the se t (|34| ) with collision kernels defined by 
formula (|l^) or j22[) and with collision rates connected 
to kernels by relation (|l5|) . The proposed proper form of 
BBE ensures that the required positivity of the atomic 
density operator is preserved. Moreover, we believe that 
the presented interpretation of the structure and terms 
of BBE will appear useful in future applications and in 
further research. 
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